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Abstract 



Wc present the analysis of the interstellar 7-ray emission measured by the Fermi Large Area 
Telescope toward a region in the second Galactic quadrant at 100° < I < 145° and —15° < b < 
+30°. This region encompasses the prominent Gould-Belt clouds of Cassiopeia, Cepheus and the 
Polaris flare, as well as atomic and molecular complexes at larger distances, like that associated 
with NGC 7538 in the Perseus arm. The good kinematic separation in velocity between the 
local, Perseus, and outer arms, and the presence of massive complexes in each of them make 
this region well suited to probe cosmic rays and the interstellar medium beyond the solar circle. 
The 7-ray emissivity spectrum of the gas in the Gould Belt is consistent with expectations based 
on the locally measured cosmic-ray spectra. The 7-ray emissivity decreases from the Gould 
Belt to the Perseiis arm, but the measured gradient is flatter than expectations for cosmic-ray 
sources peaking in the inner Galaxy as suggested by pulsars. The Xco = N{E.2)/Wco conversion 
factor is found to increase from (0.87 ± 0.05) x 10^° cm"^ (K km s"-^)"-^ in the Gould Belt to 
(1.9 ± 0.2) X 10^° cm~^ (K km s""'^)"'^ in the Perseus arm. We derive masses for the molecular 
clouds under study. Dark gas, not properly traced by radio and microwave surveys, is detected 
in the Gould Belt through a correlated excess of dust and 7-ray emission: its mass amounts to 
~ 50% of the CO-traced mass. 

Subject headings: cosmic rays - diffuse radiation - gamma rays: observations - ISM: clouds 
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1. Introduction 

Galactic interstellar 7-ray emission is produced 
through the interactions of high-energy cosmic 
rays (CRs) with the gas in the interstellar medium 
(ISM) (via pion production and Bremsstrahlung) 
and with the interstellar radiation field (via In- 
verse Compton, IC, scattering). Thus, since early 
studies with the COS-B satellite, diflPusc 7 rays 
were recognized to be a tracer of the CR densities 
and of ISM column densities in the Galaxy (Le- 
brun et al. 1983; Strong et al. 1988; Bloemen 

1989) . 

The interpretation of the observed emission is 
often based on two radio tracers of the interstellar 
gas: the 21 cm line of the hyperfine transition of 
atomic hydrogen (H i) is used to derive its column 
density A/^(Hi); the 2.6 mm line of the rotational 
transition J = 1 — ^ of CO is used to trace the 
molecular gas. The molecular phase of the ISM 
is composed mainly of H2 which cannot be traced 
directly in its cold phase. It has long been ver- 
ified, primarily using virial mass estimates, that 
the brightness temperature of CO integrated over 
velocity, Wco , roughly scales with the total molec- 
ular mass in the emitting region (see e.g. Solomon 
& Barret 1991). The conversion factor that trans- 
forms Wco into H2 column density is known as 
Xco = N(H.2)/Wco (Lebrun et al. 1983). 

The Xco conversion factor has often been as- 
sumed to be uniform across the Galaxy. We now 
have evidence, however, that it should increase in 
the outer Galaxy: from virial masses (Digel et al. 

1990) , from COBE/DIRBE studies (Sodroski et 
al. 1995, 1997) and from the measurement of the 
Galactic metallicity gradient (Israel 1997, 2000). 
A precise estimate of the Xqo gradient is neces- 
sary to measure the masses of distant II2 clouds, 
but it also impacts the derivation of the distribu- 
tion of cosmic-ray sources from 7-ray observations 
(Strong et al. 2004b). 

For many years supernova remnants (SNRs) 
have been considered the best candidates as CR 
sources. We have recently detected possible sig- 
natures of hadron acceleration in SNRs thanks to 
7-ray observations in the TeV (Albert et al. 2007; 
Aharonian et al. 2008; Acciari et al 2009) and 
GeV domain (Abdo et al. 2009d). However, the 
origin of Galactic cosmic rays is still mysterious 
and, on the other hand, the distribution of SNRs 



in the Galaxy is very poorly determined (Case & 
Bhattacharya 1998), leading to large uncertain- 
ties in the models of diffuse 7-ray emission. The 
7-ray emissivity gradient of the diffuse H i gas can 
provide useful constraints on the CR density dis- 
tribution. 

Since the Doppler shift of the radio lines allows 
kinematic separation of different structures along 
a line of sight, it is possible to constrain the 7- 
ray emissivities and the subsequent Xco ratios in 
specific Galactic regions. The performance of the 
previous 7-ray telescopes did not allow very pre- 
cise measurements beyond the solar circle (Digel 
et al. 1996, 2001). The situation has recently 
been improved with the successful launch of the 
Fermi Gamma-ray Space Telescope on 2008 June 
11. The Large Area Telescope (LAT) on board the 
Fermi mission (Atwood et al. 2009) has a sen- 
sitivity more than an order of magnitude greater 
than the previous instrument EGRET on board 
the Compton Gamma-Ray Observatory and a su- 
perior angular resolution. 

We present here the analysis of the interstel- 
lar 7 radiation measured by the Fermi LAT in 
a selected region of the second Galactic quad- 
rant, at 100° < / < 145°, -15° < 6 < -h30°, 
during the first 11 months of the science phase 
of the mission. The region was chosen because 
here the velocity gradient with Galactocentric dis- 
tance is very steep, resulting in good kinematic 
separation which allows four different regions to 
be defined along each line of sight: the nearby 
Gould Belt, the main part of the local arm, and 
the more distant Perseus and outer spiral arms. 
Among the most conspicuous clouds, one finds 
Cassiopeia, the Ccpheus and Polaris flares in the 
Gould Belt (Pcrrot & Grenier 2003; Heithausen 
&Thaddeus 1990: Grcnicr et al. 1989), the most 
massive molecular complex in the Pcrscnis arm as- 
sociated with NGC 7538 and Cas A (Ungerechts 
et al. 2000), and the off- plane molecular cloud in 
the Perseus arm associated with NGC 281 (Sato et 
al. 2007). These prominent cloud complexes are 
well suited to probe CRs and the ISM. The moti- 
vations of this work are both to provide improved 
constraints for diffuse emission models to be used 
in the detection and analysis of LAT sources and 
to reach a better comprehension of the physical 
phenomena related with diffuse 7-ray emission in 
the outer Galaxy. 
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2. Interstellar gas 

Here we describe the preparation of the maps 
tracing the column densities of the different com- 
ponents of the ISM, used in the following section 
to analyse LAT data. 

2.1. Radio and Microwave data 
2.1.1. Hi 

Column densities Af(Hi) of atomic hydrogen 
have been derived using the LAB Hi survey by 
Kalberla et al. (2005). The LSR velocity ^ cov- 
erage spans from —450 km to +400 km 
with a resolution of 1.3 km s~^. The survey angu- 
lar resolution is about 0.6°. Owing to the strong 
absorption against the radio continuum emission 
of the Cas A supernova remnant, the Hi column 
densities within 0.5° from its position were deter- 
mined by linear interpolation of the adjacent lines. 

The column densities have been derived apply- 
ing an optical depth correction for a uniform spin 
temperature Ts = 125 K, in order to directly com- 
pare our results with previous studies (like Digel et 
al. 1996). There is not general agreement in the 
literature about the values of the spin temperature 
in the atomic phase of the ISM. Prom observations 
of the 21 cm line of H i seen in absorption, Mohan 
et al. (2004a, b) derived for our region values of Tg 
varying from ~ 50 K to > 2000 K, with a mean 
value ~ 125 K. Recently Dickey et al. (2009), 
on the basis of other Hi absorption surveys, re- 
ported a mean value in the second Galactic quad- 
rant Ts = 250 K, almost constant with Galacto- 
ccntric radius. The maximum difference between 
the values of A^(H i) obtained with Ts = 125 K and 
those obtained in the optically-thin approximation 
(corresponding to the lower possible amount of gas 
or to the limit of very high spin temperature) is 
30% — 40%, whereas the maximum difference be- 
tween Ts = 250 K and optically-thin approxima- 
tion is 10% — 15%. The optical depth correction is 
non-linear, so assessing the effects of the approx- 
imation is not trivial: in particular we note that 
the uncertainties are larger where the gas density 
is higher and that assuming lower values for Ts 
we obtain structured excesses in modeled diffuse 
7-ray intensities following the shape of the clouds. 



^Local Standard of Rest velocity, i.e. the velocity in a ref- 
erence frame following the motion of the solar system. 



The systematic errors arc even larger in the 
Galactic plane where self absorption phenom- 
ena become important, especially in the Perseus 
arm where the subsequent uncertainties of derived 
iV(Hi) can reach 30% (Gibson et al. 2005). 

2.1.2. CO 

Intensities Wco of the 2.6 mm line of CO have 
been derived from the composite survey of Dame 
et al. (2001), with sampling every 0.125° near the 
Galactic plane and in the Gould Belt clouds, sup- 
plemented with observations at 0.25° sampling for 
high-latitude clouds (> 5°), covering in particular 
the region of NGC 281. 

Lines of sight not surveyed in CO were re- 
stored by linear interpolation of adjacent direc- 
tions where possible; otherwise they were assumed 
to be free of significant CO emission. CO data 
have been filtered with the moment-masking tech- 
nique in order to reduce the noise while keeping 
the resolution of the original data and retaining 
the edges of the CO clouds (see e.g. Dame et al. 
2001, Section 2.3). Preserving the faint CO edges 
is important to help decrease the degree of spa- 
tial correlation that naturally exists between the 
Af(Hi) and VFco maps of a given cloud complex 
because of the ISM multi-phase structure. 

2.2. Kinematic separation of the Galactic 

structures 

Our aim is to separately determine the 7-ray 
emission from the different Galactic structures 
present along the line of sight in the second quad- 
rant: 

1. the very nearby complexes in the Gould Belt, 
within ^ 300 pc from the solar system; 

2. the main part of the local arm, typically 1 
kpc away; 

3. the Perseus arm, 2.5 to 4 kpc away; 

4. the outer arm and beyond. 

The separation between the Gould-Belt and 
local-arm components is important to probe for 
a possible change in cosmic-ray densities between 
the quiescent nearby clouds of Cassiopeia and 
Cepheus, that produce few low-mass stars, and the 
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more active regions of the local arm which shelter 
several OB associations (Grenier et al. 1989). 

The good kinematic separation of the interstel- 
lar gas in this part of the sky is illustrated in Fig. 1. 

Galactic longitude (deg) 



145 140 135 130 125 120 115 110 105 100 




Fig. 1. — Hi longitude- velocity diagram ob- 
tained by integrating the brightness temperature 
in Kalberla et al. (2005) for |6| < 10°. The 
color scale is logarithmic in units of deg K. The 
three curves bound the preliminary Galactocentric 
rings used for analysis. At R — 8.8 kpc, i? = 10 
kpc, and R — 14 kpc (from top to bottom) they 
roughly separate Gould Belt, local arm, Perseus 
and outer arm. The separation between Gould 
Belt and main part of the local arm is hard to 
distinguish in this diagram. 

The separation of the structures along the line 
of sight was achieved through a 3-step procedure: 

a) preliminary separation based on Galactocen- 
tric rings; 

b) transformation of the ring-velocity bound- 
aries into "physical" boundaries based on 
the {I, b, v) coherence of clouds, and produc- 
tion of A^(H i) and Wco maps; 

c) correction of the A^(Hi) maps for the spill- 
over between adjacent regions. 

The three steps are described in detail below. In 
Fig. 2 the procedure is illustrated for an example 
direction at I = 133° 6 = 0°. 

The preparation of the gas maps started from 
preliminary velocity boundaries given in terms of 
Galactocentric rings that roughly encompass the 
Gould Belt for R < 8.8 kpc, the main part of the 



local arm at 8.8 kpc < i? < 10 kpc, the Perseus 
arm at 10 kpc < i? < 14 kpc, and the outer arm 
for i? > 14 kpc. Following lAU recommendations, 
we adopted a flat rotation curve with Rq — 8.5 
kpc and a rotation velocity of 220 km s~^ at the 
solar circle. The confusion that is apparent in the 
longitude- velocity {l,v) diagram of Fig. 1 between 
the Gould-Belt and local-arm components results 
from the integration over latitude and is much re- 
duced in the actual (/,&,«) cube which is used to 
construct the maps. The presence of two different 
components is evident in the example direction of 
Fig. 2: the first component peaks at ~ km s"^ 
(Gould Belt), the second one at ?; ~ —15 km s~^ 
(local arm). 

Starting from this preliminary separation, the 
ring-velocity boundaries were adjusted for each 
line of sight to better separate structures on the 
basis of their coherence in the (/, b, v) phase space. 
For each line of sight, every boundary was moved 
to the nearest minimum in the Hi spectrum, or, 
if a minimum was not found, to the nearest sad- 
dle. The shifts are typically of the order of 1 — 10 
km s~^ (see Fig. 2 b). The adjusted boundaries 
were used to calculate A^(Hi) and Wco in each 
region. 

The broad Hi clouds can easily spill over from 
one velocity interval into the next. To correct for 
this cross-contamination between adjacent inter- 
vals, for each line of sight the Hi spectrum has 
been fitted by a combination of Gaussians (see 
Fig. 2 c). The overlap estimated from the fit was 
used to correct the column density A^(Hi) calcu- 
lated in a specific interval from the spill-over from 
the adjacent regions. The correction on A^(Hi) is 
typically of the order of 1% — 10%, although it can 
reach 20% — 30% in regions corresponding to the 
frontier between clouds in the Gould Belt and in 
the main part of the local arm. 

This separation scheme provides more accurate 
estimates of the actual gas mass in a specific region 
and helps with separating structures. The result- 
ing maps are shown in Figure 3. They exhibit a 
low level of spatial degeneracy between the cloud 
complexes found in the four regions along these 
directions. Hence, we can model the observed 7- 
ray flux as a combination of contributions com- 
ing from CR interactions in the different regions. 
The correlation between the Hi and CO phases in 
each region is unavoidable, but not tight enough to 
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HI(c) 




LSR velocity (m/s) 



Fig. 2. — An example of the separation procedure described in 2.2 for the direction I — 133° 6 = 0°. Each pfot 
shows the brightness temperature Tb for H i or CO as a function of LSR velocity. Vertical lines correspond 
to the boundaries Gould Belt - local arm, local arm - Perseus arm, Perseus arm - outer arm (from right to 
left). The three rows correspond to: a) preliminary ring-boundaries, b) "physical" boundaries, c) Gaussian 
fitting of the H i line. 
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Fig. 3. — Maps of Af(Hi) (units of 10^" atoms cm~^) and Wqo (units of K km s~^). Regions: (1) Gould 
Belt, (2) local arm, (3) Perseus arm, and (4) outer arm. The maps have been smoothed for display with a 
Gaussian with a = 1°. Data sources are described in the text. 
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hamper the separation between the 7-ray emission 
from the two phases. No significant CO emission is 
found in the outer-arm region, so the correspond- 
ing map was removed from the analysis. 

2.3. Interstellar reddening 

An excess of 7 rays (observed by EGRET) 
correlated with an excess of dust thermal emis- 
sion was found over the 7V(Hi) and Wco column- 
density maps in all the nearby Gould-Belt clouds 
by Grenier et al. (2005). Therefore, they reported 
a considerable amount of "dark" gas, i.e. neutral 
gas not properly traced by Hi and CO, at the in- 
terface between the two radio-traced phases. The 
chemical state of the additional gas has not been 
determined yet, leaving room for H2 poorly mixed 
with CO or to H i, overlooked e.g. because of in- 
correct assumptions about the spin temperature 
for optical depth corrections or H i self absorption 
(see 2.1.1). 

Following the method proposed by Grenier et 
al. (2005), we have prepared a map to account 
for the additional gas. The map is derived from 
the E(B — V) reddening map of Schlegel et al. 
(1998), which provides an estimate of the total 
dust column densities across the sky. Point sources 
(corresponding to IRAS point sources) were re- 
moved and the corresponding pixels were set to 
the average value of the adjacent directions. In 
order to subtract the dust components correlated 
with iV(Hi) and Wqo, the reddening map was 
fitted with a linear combination of the same set 
of A''(Hi) and Wco maps for the Gould-Belt and 
local, Perseus, and outer-arm regions described 
above. A detailed discussion of the results of the 
fit goes beyond the scope of the present work, 
so it is deferred to another paper (A. A. Abdo 
et al. 2009, in preparation), which will address 
the results over several interstellar complexes in 
the Gould Belt and will compare them with 7-ray 
measurements by the LAT. 

The resulting E(B — V)^.^,^ residual map, ob- 
tained subtracting from the E(B — V) map the 
best-fit linear combination of our set of iV(H i) and 
Wco maps, is shown in Fig. 4. 

The residuals typically range from —1 to +1 
magnitude. Unlike in Grenier et al. (2005), both 
positive and negative residuals have been consid- 
ered in the analysis of the 7-ray data. Residu- 
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Fig. 4. — E(B — V)j.(,g map: map of the reddening 
residuals obtained after subtraction of the parts 
linearly correlated with the combination of Af(Hi) 
column densities and Wco intensities found in the 
four regions along the line of sight (Gould Belt, lo- 
cal, Perseus, and outer arms). The positive resid- 
uals surrounding CO clouds off the plane outline 
the potential dark gas envelopes of the Gould Belt 
clouds. The map has been smoothed for display 
with a Gaussian with cr = 1°. 

als hint at limitations in the gas radio tracers as 
well as in the E(B — V) map. Positive residuals 
can correspond to a local increase in the dust- 
to-gas ratio and/or to the presence of additional 
gas not properly accounted for in the Af(Hi) and 
Wco maps. The latter explanation is supported 
by the significant correlation we will find between 
the E(B — V)j.j,g map and the LAT 7-ray data 
(see 4.2.2). Fig. 4 shows that at |6| > 5° the 
E(B — V)j.j,g map is dominated by positive resid- 
uals forming structured envelopes around the CO 
clouds. Small negative residuals are systematically 
seen toward the CO cores. They may be due to 
a decrease in dust temperature in the denser, well 
shielded, parts of the molecular complexes, or to 
local variations of the dust-to-gas ratio. Further 
comparison with dust extinction tracers is needed 
to investigate this effect. Positive and negative 
residuals appear at low latitude, but, because of 
the pile up of dust clouds with different tempera- 
tures along the line of sight, the temperature cor- 
rection, applied by Schlegel et al. (1998) to the 
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thermal emission to produce the E(B — V) map, is 
not as valid near the plane as in well resolved lo- 
cal clouds off the plane, in particular toward bright 
star-forming regions. These effects may cause the 
clusters of negative residuals at \b\ < 5°. We note 
that the positive residuals along the plane are not 
well correlated with the amount of self absorbed 
Hi found in the second quadrant (Gibson et al. 
2005). The most conspicuous self-absorbed Hi 
cloud in fact corresponds to the negative residuals 
seen at 130° < I < 140°. 

By construction, the E(B — V)^j,g map compen- 
sates for the limitations of the radio tracers, both 
by revealing non-emitting gas and by correcting 
the approximations applied to handle the radia- 
tive transfer of the radio lines. In particular the 
E(B — V)j.j,g map depends on the optical depth 
correction applied to the A^(Hi) maps. We note 
that (lifFcrcnt choices of the H I spin temperature, 
from the optically thin case down to 250 and 125 
K, change the E(B — V) residuals by up to 0.2—0.3 
mag on the plane and 0.1 mag at |fe| > 5°. Off the 
plane, the shape of the large structures of posi- 
tive residuals around the Gould-Belt clouds hardly 
changes. 

The strong correlation between E(B — V) resid- 
uals and 7-ray data, which will be shown in 4.2.2, 
proves that the interstellar reddening is in many 
directions a better tracer of the total gas column 
densities than the combination of Hi and CO. 
Therefore, we will use the E(B — V)^.^^ map to 
correct the standard radio and microwave tracers, 
very suitable for the aims of this work since, unlike 
reddening, they carry distance information. 

3. Gamma-ray analysis 
3.1. LAT data 

The Fermi LAT is a pair-tracking telescope 
(Atwood et al. 2009), detecting photons from 
- 20 MeV to more than 300 GeV. The tracker 
has 18 {x, y) layers of silicon microstrip detectors 
interleaved with tungsten foils to promote the con- 
version of 7 rays into electron-positron pairs (12 
thin foils of 0.03 radiation lengths in the front sec- 
tion plus 4 thick foils of 0.18 radiation lengths in 
the back section; the last two layers have no con- 
version foils). The tracker is followed by a seg- 
mented Csl calorimeter to determine the 7-ray en- 
ergy. The whole system is surrounded by a scintil- 



lator shield to discriminate the charged cosmic-ray 
background. The instrument design and the anal- 
ysis result in a peak effective area of ~ 8000 cm^ 
(~ 6 times greater than EGRET), a field of view 
of ~ 2.4 sr (~ 5 times greater than EGRET) and a 
superior single photon angular resolution (for front 
converting photons, the 68% containment angle at 
1 GeV reaches ~ 0.6° with respect to ~ 1.7° for 
EGRET). 

Data were obtained during the period 2008 Au- 
gust 4 - 2009 July 4. The Fermi observatory was 
operated in scanning sky survey mode, rocking 
35° north and south of the zenith on alternate or- 
bits, apart from calibration runs that are excluded 
from the analysis. We used the dataset prepared 
for the construction of the first year Catalog of 
LAT sources (A. A. Abdo et al. 2009, in prepara- 
tion), excluding brief time intervals corresponding 
to bright 7-ray bursts. It uses the Diffuse event 
selection, which has the least residual CR back- 
ground contamination (Atwood et al. 2009). We 
also selected events on the basis of the measured 
zenith angle to limit the contamination from inter- 
actions of cosmic rays with the upper atmosphere 
of the Earth. Owing to these interactions the limb 
of the Earth is a very bright 7-ray source, seen at 
a zenith angle of ~ 113° at the 565 km nearly- 
circular orbit of Fermi. Since our region is close 
to the North celestial pole it is often observed at 
large rocking angles. In order to reduce the Earth 
albedo contamination, we accept for analysis here 
only events seen at a zenith angle < 100°. The 
exposure is only marginally affected (because the 
detection efficiency dramatically decreases at large 
inclination angles), but the background rate is sig- 
nificantly reduced. 

3.2. Model for analysis 

The analysis scheme used since the COS-B era 
(Lebrun et al. 1983; Strong et al. 1988; Digel 
et al. 1996) is based on a very simple transport 
model. Assuming that the interstellar medium is 
transparent to 7 rays, that the characteristic diffu- 
sion lengths for CR electrons and protons exceed 
the dimensions of cloud complexes, and that cos- 
mic rays penetrate clouds uniformly to their cores, 
the 7-ray intensity / (cm~^ st~^) in a direction 
{I, b) can be modeled to first order as a linear com- 
bination of contributions coming from CR interac- 
tions with the different gas phases in the various 
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regions along the line of sight. Wc add the con- 
tribution from point-like sources and an isotropic 
intensity term. Several processes are expected to 
contribute to the latter, notably the cxtragalac- 
tic 7-ray background and the residual instrumen- 
tal background from misclassified interactions of 
charged CRs in the LAT. The IC emission is also 
expected to be rather uniform across this small re- 
gion of the sky. We used the current best models 
of IC emission to verify that it is statistically not 
distinguishable from an isotropic background over 
the small region of interest, at large angular dis- 
tance from the inner Galaxy (sec 3.4). The present 
analysis does not aim to provide meaningful re- 
sults for the extragalactic background and the IC 
emission which will be addressed in forthcoming 
publications (Abdo et al. 2009f, A. A. Abdo et 
al. 2009, in preparation). 

In the absence of suitable tracers for the dif- 
fuse ionized gas (primarily Hii), the derived 7-ray 
emissivities for neutral gas will be slightly over- 
estimated. However, the ionized gas is contribut- 
ing to ~ 10% of the total mass and, because of 
its large scale height of ~ 1 kpc above the plane 
(Cordes & Lazio 2002), part of its 7-ray emission 
will be overtaken by the isotropic term in the fit 
to the LAT data. So the bias on the neutral gas 
emissivities should be small. 

Therefore, the 7-ray intensity /, integrated in a 
given energy band, is modeled by Eq. 1. 

4 

I{l,b) = ^ [(?H M • A^(H I) (Z, 6), + gco,»-l^co + 

2=1 

+ 9EBV • E(B - V)^^,g(Z, b) + 

+ ^5,-(5(2)(;-;^,6-6,) (1) 
J 

The sum over i represents the combination of the 
four Galactic regions. The free parameters are the 
emissivities of Hi gas, gHi,j (s~^ sr~^), per unit 
of Wco intensity, qco,i (cm~^ s~^ sr~^ (K km 
s""'^)""'^), and per unit of E(B — V) residuals, ^ebv 
(cm~^ s~^ sr^^ mag~^). /iso (cm~^ sr^^) is 
the isotropic background intensity. The contribu- 
tion from point sources is represented by the sum 
over J, where Sj is the integrated flux (cm~^ s~^) 
of the source lying at the position {Ij, bj). 



3.3. Analysis procedure 

3.3.1. Method 

We used the standard LAT analysis environ- 
ment provided by the Science Tools'^. The 7- 
ray statistics are large enough to model the spec- 
tral shape of each component as a power law in 
relatively narrow energy bands. This assump- 
tion, together with the iterative procedure de- 
scribed below in 3.3.2, allows the exposures and 
the convolution with the energy-dependent Point 
Spread Function (PSF) to be computed without 
forcing an a priori spectral index. The Sciejice 
Tools provide a full convolution of the maps with 
the energy-dependent PSF. The Science Tools 
are also very flexible in the description of point 
sources (number, location, spectra). We used the 
P6_V3 post launch Instrument Response Functions 
(IRFs), which take into account the loss of detec- 
tion efficiency due to pile-up and accidental coin- 
cidence effects in the LAT (Rando et al. 2009). 

LAT data have been analysed using a binned 
maximum-likelihood procedure with Poisson statis- 
tics, on a spatial grid with 0.5° spacing in Carte- 
sian projection. The higher energy range we have 
investigated starts at a few GeV, where the 68% 
containment angle is ~ 0.5° for events converting 
in the front section of the tracker (about a factor 
two larger for back converting events), so we can- 
not resolve details smaller than this in the 7-ray 
maps. This resolution is commensurate with that 
of the Hi and E(B — V) maps. 

The analysis was performed for 5 contiguous 
energy bands: 200 MeV - 400 MeV, 400 MeV - 
600 MeV, 600 MeV - 1 CeV, 1 GeV - 2 GeV and 
2 GeV - 10 GeV. The energy bands were chosen 
wide enough to obtain stable results for the fit 
parameters, because large statistical fluctuations 
might hamper the separation of the different maps. 
Below 200 MeV the broad PSF does not allow an 
effective separation of the different maps. We are 
confident that between 0.2 and 10 GeV the inter- 
stellar 7-ray emission from the gas dominates over 
the instrumental foregrounds. The count maps in 
the five energy bands are shown in Fig. 5. 



^ http: / / fermi.gsfc.nasa.gov/ssc / data/ analysis / documentation/Cicerone / 
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3.3.2. Point sources 

The inclusion of sources in the analysis model 
is a non-trivial task because the likelihood maxi- 
mization procedure (based on the optimization en- 
gine Minuit2^) is stable up to a few tens of free 
parameters. The sources have thus been added 
following an iterative procedure. 

The sources were taken from the 11 month 
source list, which will be the basis for the First 
Year Catalog of LAT sources in preparation 
(A. A. Abdo et al. 2009, in preparation). The 
sources were added following the detection signif- 
icance (TS) in the 11 month source list''. The 
sources are added as point-like sources keeping 
their positions at those given in the list while 
letting their power-law spectra to vary indepen- 
dently in each energy band. No further attempts 
are made in this analysis to improve the spectral 
modeling or to account for possible extension. 

The inclusion of the sources went through the 
following steps, where the parameters of the dif- 
fuse emission model were always let free: 

• We started with no point sources in the 
model. 

• We added 9 sources detected with TS > 600 
(hereafter bright sources). They were added 
to the sky model 3 at a time in order of de- 
creasing TS, freezing at each step the previ- 
ous source spectra and fitting the last three, 
while the diffuse parameters were always let 
free. Among bright sources, for the 6 sources 
lying in the region under study, we let their 
fluxes and spectral indexes free; for 3 soTirces 
lying just outside (< 5°) the region bound- 
aries, we fixed their parameters at the val- 
ues determined in the 11 month source list. 
These bright sources were already reported 
in the LAT Bright Source List (Abdo et 
al. 2009b): two of them are firmly iden- 
tified as pulsars (OFGL J0007.4 + 7303 or 
LAT PSR J0007-h7303, and OFGL J2229.0-h 
6114 or PSR J2229+6114), one as a 7-ray bi- 
nary (OFGL J0240.3-F6113or LSI-K61 303), 
and the others are associated with blazars. 



^http:/ /wwwasdoc. web. cern.ch/wwwasdoc/minuit/minmain. html 

"The Test Statistic, TS, is defined as TS = 2(liiif - 
ln.i?o), where .if and ^0 are the maximum-UkeUhood val- 
ues reached with and without the source, respectively. 



• We then added 52 more sources in the 11 
month source list within the region bound- 
aries with TS between 600 and ~ 25 (out 
of them 22 where detected with TS > 100); 
they were added in several groups of 6 or 5 
sources, with a procedure analogous to that 
used to handle bright sources, but only their 
integrated fluxes were allowed to be free, 
whereas the spectral indexes were fixed at 
the values in the source list. 

• Finally, the analysis was repeated with all 
the sources, letting free only the parame- 
ters of the diffuse model and of the bright 
sources. 

The iterative procedure allows verifying that 
only the bright sources can affect the diffuse pa- 
rameters: the latter do not significantly change 
when less significant sources (TS < 600) are added 
to the model. This docs not apply to /iso and 
In 1,4- we note that their values keep decreasing 
as we add new sources down to TS 25. We ar- 
gue that the isotropic intensity generally absorbs 
point sources off the plane that are not included 
in the analysis; as said before we are not trying to 
give a physical interpretation of /iso- On the other 
hand, given the low linear resolution in distant 
clouds of the outer arm and the subsequent lack of 
pronotmced features in the map (see Fig. 3), point 
sources at very low latitude (|6| < 3°) can strongly 
bias the value of the corresponding H i emissivity, 
qm.i, as separating them from the clumpy ISM 
emission near the plane in 0.5° maps is difficult. 
Therefore, we consider this parameter only as an 
upper limit to the real gas emissivity in the outer 
arm. 

3.4. Fit results 

The quality of the final fits is illustrated in the 
residual maps of Fig. 6. The residuals, i.e. ob- 
served counts minus model-predicted counts, are 
expressed in standard deviation units (square root 
of model-predicted counts). The maps show no ex- 
cesses below — 4(T or above +7a. 

The best-fit parameters obtained in the five en- 
ergy bands are given in Table 1, where the uncer- 
tainties correspond only to statistical errors. We 
have also evaluated the systematic errors due to 
the uncertainties on the event selection efficiency. 
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Fig. 6. — Gamma-ray residual maps in the same energy bands as in Fig. 5. The residuals, i.e. observed 
counts minus model-predicted counts, are in units of the square root of the model-predicted counts (truncated 
between —4 and +4 for display). 
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Table 1 

Parameters of the diffuse emission model obtained from the fit to LAT data. 



para.motcr''' 


' 0.2 0.4 GcV 


0.4 0.(i GoV 


0.() 1 GcV 


1 2 GcV 


2 10 Gc-V 


9Hi,1 


0.584±0.011 


0.224±0.008 


0.168±0.004 


0.110±0.003 


0.048±0.002 


QCO, 1 


1.09±0.04 


0.367±0.017 


0.318±0.013 


0.198±0.008 


0.102±0.005 


9Hi, 2 


0.536±0.018 


0.200±0.007 


0.157±0.005 


101±0.004 


0.054±0.002 


qco,2 


1.67±0.17 


0.47±0.06 


0.44±0.04 


0.26±0.03 


0.087±0.014 


9Hi,3 


0.349±0.011 


0.128±0.004 


0.108±0.003 


0.072±0.002 


0.0397±0.0014 


qco,3 


1.17±0.15 


0.52±0.06 


0.37±0.04 


0.24±0.03 


0.115±0.016 


Qui, 4: 


0.33±0.04 


0.101±0.017 


0.114±0.013 


0.103±0.009 


0.032±0.005 


Qebv 


16.7±1.0 


6.0±0.4 


3.49±0.27 


2.28±0.18 


0.80±0.11 




4.67±0.10 


1.19±0.04 


0.92±0.03 


0.63±0.02 


0.371±0.0017 



'^Units: quj,t (lO'^^ s'^ sr"!), qco,^ {W^ cm-^ s'^ sr"! (K km s-^)"!), gEBV (10"^ cm-^ 
sr~^ mag"-^), liso (10~^ cm~^ s~-^ sr~^). 

'^The subscripts refer to the different regions under analysis: 1) Gould Belt, 2) local arm, 3) 
Perseus arm, 4) outer arm and beyond. 



From the comparison between Monte Carlo sim- 
ulations and real observations of the Vela pulsar, 
they are evaluated to be 10% at 100 MeV, 5% 
at 500 MeV and 20% at 10 GcV, scaling linearly 
with the logarithm of energy between these val- 
ues. These uncertainties were parametrized into 
two sets of IRFs encompassing the most extreme 
scenarios. The last step of the analysis has been 
repeated using these two IRF sets and the results 
are assumed to bracket the systematic errors due 
to the event selection efficiency (shown as shaded 
grey areas in the following figures). Only the last 
step was considered, because we previously veri- 
fied that only bright sources impact the parame- 
ters of the diffuse emission model. 

We also verified the impact of the isotropic ap- 
proximation for the IC emission, repeating the last 
step of the analysis including a recent model based 
on the GALPROP CR propagation code (see e.g. 
Strong & Moskalenko 1998; Strong et al. 2004a; 
Porter et al. 2008). The values obtained for 
the parameters of the diflfuse emission model were 
compatible with the previous results, except for 
the isotropic intensity. 

Other systematic uncertainties will be ad- 
dressed in the discussion section. 



4. Discussion 

4.1. Emissivity per Hi atom and cosmic- 
ray spectra 

4.1.1. Consistency with other measurements 

In Fig. 7 we report the emissivity spectra per 
Hi atom measured in the Gould Belt, the main 
part of the local arm and the Perseus arm. The 
inclusion of the E(B — V)j.gg map in the fit does 
not have a strong impact on the emissivities of the 
broadly distributed H I gas, which decrease by less 
than 10% considering the interstellar reddening in 
analysis. 

The results we obtained in the Gould-Belt and 

local-arm regions are consistent below 1 GeV with 
the measurement by Digel et al. (1996), ob- 
tained from EGRET observations of the region of 
Cepheus and Polaris. Above 1 GeV LAT measure- 
ments are ~ 40% lower than those by EGRET. We 
will see below in 4.1.3 that LAT measurements 
are consistent with the a priori expectations for 
the local Hi emissivity: this result confirms that, 
as was already deduced from LAT observations of 
broader regions of the sky (Porter et al. 2009; 
Abdo et al. 2009e), LAT measurements are not 
consistent with the GeV excess seen by EGRET, 
which was noticed also as an excess above 1 GeV 
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Fig. 7. — Emissivity spectra per Hi atom as measured in the Gould Belt (1), local arm (2), and Perseus 
arm (3) clouds. Horizontal bars mark the energy bands, vertical bars show the statistical uncertainties on 
the measurement. The shaded areas represent the systematic errors due to the uncertainties on the event 
selection efficiency. 



in the emissivity of nearby Hi complexes, as dis- 
cussed in Digel et al. (2001). 

Our spectra of the emissivity per Hi atom are 
consistent with the results of an independent anal- 
ysis carried out on LAT data to determine the lo- 
cal Hi emissivity in a mid-latitude region of the 
third Galactic quadrant (Abdo et al. 2009c). The 
latter analysis investigated a different region of the 
sky, but encompassing H i complexes at < 1 kpc 
from the solar system, mostly located in a seg- 
ment of the local arm. If we compare the present 
results in the main part of the local arm with those 
by Abdo et al. (2009c) we have excellent agree- 
ment. Therefore we have verified that CR proton 
densities smoothly vary on a few kpc scale around 
the solar system. 

4-. 1.2. Physical model 

We further compare our results with the pre- 
dictions by GALPROP, a physical model of CR 



propagation in the Galaxy (see e.g. Strong & 
Moskalenko 1998; Strong et al. 2004a, 2007). 
GALPROP solves the propagation equation for all 
CR species, given a CR source distribution and 
boundary conditions. Current GALPROP mod- 
els assume a Galactocentric source distribution de- 
rived from that of pulsars (Strong et al. 2004b). 
The distribution used by the model adopted for 
this work, called 54_71Xvarh7S, is given by Eq. 2 



exp 



R ^ R, 



Rr. 



(2) 



with a = 1.25, f3 = 3.56 and Rq = 8.5 kpc. A 
truncation is applied at i? = 15 kpc because we 
do not expect many CR sources in the outermost 
Galaxy. This choice of parameters results in a 
slightly flatter radial profile of CR densities than 
with the pulsar distribution. 

The GALPROP model 54_71Xvarh7S is tuned 
to reproduce the in situ measurements of CR spec- 
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tra at the solar circle. The proton spectrum is de- 
rived from a compilation of direct measurements 
(Alcaraz et al. 2000; Sanuki et al. 2000; Shikaze 
et al. 2007). The model includes the CR electron 
spectrum recently measured by the LAT (Abdo et 
al. 2009a). 

Once the propagation equation is solved, GAL- 
PROP computes the emissivity for stable secon- 
daries, in particular 7 rays. The electron Bremm- 
strahlung component is evaluated using the for- 
malism by Koch & Motz (1959) as explained in 
Strong et al. (2000). The emissivity due to p-p 
interactions is evaluated using the inclusive cross 
sections as parametrized by Kamae et al. (2006). 
Following the method by Dermer (1986a, b), the p- 
p emissivity is increased to account for interactions 
involving CR a particles and interstellar He nuclei. 
This method provides an effective enhancement 
with respect to pure p-p emissivities, often named 
the nuclear enhancement factor, of en — 1.45. 
More recent calculations by Mori (2009), how- 
ever, report values as large as en — 1.75 — 2 due 
to different CR spectral formulae (Honda et al. 
2004) , different ISM abundances and the inclusion 
of heavier nuclei both in CRs and in the ISM. Fur- 
ther theoretical developments are required to bet- 
ter constrain en, extending the predictions from 7 
rays to other relevant messengers like antiprotons 
(Adriani et al. 2009). 

4.. 1.3. Emissivity in the Gould Belt 

In Fig. 8, we compare the emissivity spectrum 
per H I atom we measured in the Gould Belt with 
the GALPROP predictions. We find the latter to 
be ~ 50% lower at all energies. As we have just 
discussed, a large part of this excess (~ 30%) can 
be explained by the uncertainties in the contribu- 
tion from interactions involving CR and ISM nu- 
clei other than protons. The remaining ~ 20% ex- 
cess can be explained by systematic uncertainties 
in the CR proton spectra at the Earth (~ 20%), 
the A^(Hi) column-density derivation, and the 
kinematical separation of emission from the outer 
Galaxy. 

The Hi emissivity in the Gould Belt clouds 
(within 300 pc from the solar system) is thus con- 
sistent with the hypothesis that the gas is interact- 
ing with CRs with the same spectra measured at 
Earth. Fig. 8 shows the GALPROP model scaled 
by -f-50% to highlight that the spectral shape is in 



■ Bremsstrahlung 

■ pion decay 

■ total 
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Fig. 8. — Hi emissivity spectrum in the Gould 
Belt, as shown in Fig. 7. The curves represent 
the predictions by GALPROP 54_71Xvarh7S. The 
total emissivity from the model has been increased 
by 50% to reproduce our measurements. 

good agreement with our results. 

The Hi spin temperature of 125 K (chosen to 
have a straightforward comparison with earlier 
analyses) is among the lowest values reported in 
the literature. A higher temperature would imply 
a higher emissivity, therefore a larger discrepancy 
with the GALPROP model (e.g. gni. 1 increases 
by another 5% to 10% if we take Ts = 250 K as 
recently suggested by Dickey et al. 2009). 

4-. 1.4- Hi emissivity gradient 

It is evident from Fig. 7 that the Hi emissiv- 
ity decreases from the Gould Belt to the Perseus 
arm, as expected from the declining distribution of 
candidate CR sources in the outer Galaxy. Fig. 9 
shows the emissivity ratios between the more dis- 
tant regions and the Gould Belt. Systematic er- 
rors due to the event selection efficiency are not 
relevant for these ratios, because the emissivity 
spectra are similar. 

The emissivity spectrum in the local arm is 
10% lower than in the Gould Belt. The GAL- 
PROP model predicts such a decrease because of 
the change in Galactocentric radius from the solar 
circle to the main part of local arm, located in this 
direction at ~ 9.5 kpc. 

A further decline in i is expected between the 
local and Perseus arms, but Fig. 9 shows that the 
measured gni gradient is significantly shallower 
than the GALPROP prediction. In Fig. 10 we 
compare the Hi emissivity integrated above 200 
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Fig. 9. — Emissivities measured in the Perseus 
arm -black- and in the local arm -light (red)- rel- 
ative to those in the Gould Belt. Horizontal bars 
mark the energy ranges, vertical bars show statis- 
tical uncertainties. The dashed lines represent the 
GALPROP predictions. 
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nation is probably affected by faint sources (see 
section 3.3.2). 

The discrepancies between the measured and 
predicted gradients may be due to the large un- 
certainty in the CR source distribution. The SNR 
radial distribution across the Galaxy is very poorly 
determined because of the small sample available 
and large selection effects (Case & Bhattacharya 
1998). Distance and interstellar dispersion uncer- 
tainties also bias the pulsar distribution, in spite 
of the larger sample available (Lorimer 2004). On 
the other hand, the CR diffusion parameters, de- 
rived from local isotopic abundances in CRs, may 
not apply to the whole Galaxy, as suggested by 
Taillet & Maurin (2003). Self absorption can also 
lead to a significant underestimate of (H i) in the 
Perseus arm (Gibson et al. 2005), thus to an over- 
estimate of its 7-ray emissivity. Therefore, further 
investigation is needed to better understand the 
radial profile of the H i emissivity. 

In Fig. 7, the H i emissivity spectrum in the 
Perseus arm appears harder than expectations, 
thus suggesting that primary CR spectra vary 
across the Galaxy. We cannot, however, rule out 
energy-dependent systematic effects due to the 
separation power provided by the LAT PSF which 
strongly varies with energy, or a hardening due to 
contamination by hard unresolved point sources, 
like pulsars, clustering in the Perseus-arm struc- 
tures. 



Fig. 10. — Radial profile with Galactocentric ra- 
dius of the H I emissivity integrated between 200 
MeV and 10 GeV. Black dots/horizontal bars 
mark the ranges in kinematic distance encompass- 
ing the Gould Belt, the main part of the local arm, 
the Perseus and outer arms (from left to right). 
Statistical uncertainties on qui are smaller than 
the dot dimensions. The grey shaded area shows 
the systematic uncertainties on the event selection 
efficiency. The (blue) dashed line is the GAL- 
PROP prediction scaled up by 50%. 

MeV predicted by GALPROP as a function of 
Galactocentric radius with the values we measured 
in the four regions defined for analysis, drawing 
the same conclusion. In this figure we report the 
emissivity found in the outer arm, though consid- 
ered only as an upper limit because its determi- 



4.2. Cloud masses 
12.1. CO 

Because the 7-ray emission from molecular 
clouds is primarily due to H2 and the molecu- 
lar binding energy is negligible with respect to 
the energy-scale of the 7 radiation processes, the 
emissivity per H2 molecule is twice the emissivity 
per H I atom. Under the hypothesis that the same 
CR fiux penetrates the Hi and CO phases of a 
cloud, we can assume that qco,i = 2 ATcci ■ ^Hi.i 
in each region to derive the CO-to-H2 conversion 
factor, Xco- 

We have performed a maximum likelihood lin- 
ear fit qco,2 = A:co,i-29hi,i+9, between the qco,i 
and qn i, i values found in the various energy bands 
for each region. We have taken into account the 
errors and covariances obtained from the 7-ray fits 
for both qn i and qco ■ Systematic errors due to the 
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Fig. 11. — Correlation between the Hi and CO emissivities obtained in the five energy bands for each region 
under analysis: (1) Gould Belt, (2) local arm, (3) Perseus arm. Error bars show the statistical uncertainties 
on (7hi,j and qco,i- The (red) lines give the best linear fits. 



event selection efficiency do not affect the deriva- 
tion of the Xco slope because the Hi and CO 
emissivities have similar spectra. The results are 
shown in Fig. 11, and the best-fit parameters are 
reported in Table 2. We observe a good linear cor- 
relation between gni and gco that lends support 
to the assumption that CRs penetrate molecular 
clouds uniformly to their cores (still under debate, 
see e.g. Gabici et al. 2007). 

Fig. 12 shows the Xco variation with Galacto- 
centric radius. Our measurements are consistent 
with previous 7-ray estimates in this region of the 
sky (Digel et al. 1996), but they are more pre- 
cise, especially in the outer Galaxy. For the seg- 
ment of the Perseus arm near NGC 7538 we have 
lowered the statistical uncertainty from ^ 40% to 
10%. The results suggest an increase of Xqo in 
the outer Galaxy, as expected from the metallicity 
gradient (see e.g. Rolleston et al. 2000). The Xqo 
measurements in external galaxies indeed show 
a metallicity dependence possibly caused by CO 



photodissociation and poor self-shielding in low- 
metallicity environments (Israel 1997, 2000). 

Contamination from unresolved point sources 
with a spatial distribution closely related to that 
of the clouds is expected in star-forming regions 
which can produce young pulsars, supernova rem- 
nants, and massive binaries. This effect is unlikely 
in the Gould-Belt clouds (Cassiopeia, Cepheus, 
and Polaris), first because they form few high- 
mass stars, second because of the good linear res- 
olution of the 7-ray maps of these nearby clouds. 
Their proximity (< 300 pc) and the ~ 0.5° angular 
resolution of the LAT in the higher energy band 
imply a linear resolution < 3 pc, which allows an 
efficient separation between diffuse emission and 
point sources. The contamination by point sources 
is limited for similar reasons in the nearby local 
arm, < 1 kpc away, but it cannot be clearly ruled 
out in the Perseus arm clouds which are known 
to form massive star clusters (see e.g. Sandell & 
Sievers 2004). 
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Table 2 

Results of the linear fits between the Hi and CO emissivities in the different regions (1 
- Gould Belt, 2 - local arm, 3 ~ Perseus arm): (jco,i = Xcoi ■ 2gHi,i + g^- 



1 0.87±0.05 0.015±0.012 

2 1.59±0.17 -0.08±0.03 

3 1.9±0.2 -0.03±0.03 



'^Units: lO^o cm-^ (K km 
s-i)-i). 

''Units: 10"*^ cm-^ s'^ sr"! 
(K km s-^)-i. 



'" 2 4 6 8 10 12 14 

R (kpc) 

Fig. 12. — Xco as a function of Galactocentric 
radius. The solid-line (black) points represent our 
measurements: horizontal bars mark the ranges of 
kinematic distance encompassing the Gould Belt 
and the local and Perseus arms (from left to right), 
vertical bars show the statistical uncertainties on 
Xco (errors are statistical only, possible system- 
atics are discussed in the text in 4.2.1). Dashed 
(blue) lines represent the values used in GAL- 
PROP by Strong et al. (2004b). The solid 
(green) line shows the Xqq function determined 
by Nakanishi & Sofue (2006) from CO data and 
virial masses (adapted to the rotation curve as- 
sumed for our analysis). The dashed (grey) points 
show the previous EGRET measurements in the 
region of Cepheus and Polaris (Digel et al. 1996). 



We cannot exclude separation problems be- 
tween the 7-ray emission from the CO cores and 
their surrounding H I envelopes. The separation. 



based on the spatial distribution of the different 
phases, becomes less efficient with increasing dis- 
tance due to the lower linear resolution. More- 
over, we have verified that the presence of 7 rays 
associated with the dark-gas envelopes around the 
CO cores affects the determination of the CO-to- 
H2 factor in more distant, not so well resolved, 
clouds (whereas the impact is negligible in the 
closer clouds). Excluding the E(B — V)^.^,^ map 
from the model yields a ~ 30% increase of Xqo in 
the Perseus arm. Unfortunately, the E(B — V)^.^,^ 
map is not reliable near the plane because of the 
confusion along the line of sight that prevents an 
adequate temperature correction and the removal 
of clusters of IR point sources, so this difference 
has to be considered as a systematic error on Xqq 
in the Perseus arm. 

Whether the present Xqo gradient can be fully 
attributed to the metallicity gradient, or partially 
to unresolved sources, H i and CO separation prob- 
lems, or gas not traced by Hi and CO, needs 
further investigation, primarily at higher resolu- 
tion when more high-energy LAT data become 
available to profit from the better angular reso- 
lution. For the moment, the fact that the present 
Xco determination does not depend on energy 
(see Fig. 11) suggests that unresolved sources and 
separation of the different gas phases do not signif- 
icantly influence the result. The results shown in 
Fig. 12 indicate significantly smaller Xco values in 
the outer Galaxy than those used by Strong et al. 
(2004b) in GALPROP and systematically smaller 
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values than the Xco{R) relation determined by 
Arimoto et al. (1996) and Nakanishi & Sofue 
(2006) using CO data and virial masses. 

The Xco values shown in Table 2 have been 
used to estimate cloud masses using Eq. 3 

M = 2 nmucf Xco J Wco {I, h) dO (3) 

where d is the distance of the cloud, mu is the 
H atom mass, and n = 1.36 is the mean atomic 
weight per H atom in the ISM. We did not use the 
kinematic distances inferred from CO surveys, but 
we adopted more precise estimates available in the 
literature. The results are given in Table 3. The 
errors include only the statistical uncertainties on 

To investigate the discrepancies found between 
the different determinations of Xco we calculated 
the virial masses for well-resolved clouds off the 
plane. The virial masses have been obtained from 
the CO velocity dispersion for a spherical mass 
distribution with density profile oc 1/r, following 
Eq. 4 

M = Ya,l (4) 

where r is the cloud radius, (t„ the velocity dis- 
persion and G is Newton's constant. The veloc- 
ity dispersion has been measured for each line of 
sight and the average value in the sample has been 
taken as the characteristic in the cloud. This 
method limits the impact of the obvious veloc- 
ity gradients in these clouds. Bcc;aiise the virial 
mass heavily depends on the estimate of the char- 
acteristic radius and on the cut-off applied in its 
evaluation, we considered both the effective radius 
ta = -^/A/n (where A is the geometrical area 
of the cloud) and the intensity weighted radius 
(r) = (Ej Wcoz n) / (E, Wcoz) (where r, is the 
distance of pixel i to the peak Wco pixel). We 
truncated the calculation at 1% of the Wco peak 
in both cases. We find that the virial masses are 
systematically larger than the Xco derived masses 
by a factor 1.5-3. This discrepancy in the nearby 
clouds is comparable to that shown in Fig. 12 be- 
tween the 7-ray estimates of Xco and the Xco (-R) 
function by Nakanishi & Sofue (2006) which relies 
on virial masses. The 7-ray estimates are indepen- 
dent from the chemical, dynamical and thermo- 
dynamical state of the clouds, but they can suf- 
fer from the limited resolution of 7-ray surveys 



and the non-uniform penetration of CRs into the 
dense CO cores. Conversely, the assumption of a 
spherical cloud in virial equilibrium against tur- 
bulent motions is rather crude. Intrinsic velocity 
gradients and magnetic pressure can easily bias 
the virial mass results. 

4-2.2. Dark gas 

In order to quantify the significance of the 
correlation between the 7-ray intensities and the 
E(B — V)j.j,g map, we have repeated the last step 
of 3.3.2 without including it in the analysis. The 
corresponding test statistics, TS = 2A(ln^), ob- 
tained in the five energy bands are given in Ta- 
ble 4. With the addition of two free parameters 
(^EBV and a spectral index), in the null hypothesis 
that there is no 7-ray emission associated with the 
E(B - V)^gg map TS should follow a distribu- 
tion with two degrees of freedom. Therefore, the 
correlation between 7 rays and E(B — V) residu- 
als is verified at a confidence level > 99.9% in all 
energy bands. 

The magnitudes of the dust masses and dust 
IR emission are too low to explain this correlation 
by CR interactions with dust grains or their ther- 
mal radiation. However, the correlation can be 
explained by CR interactions in normal gas that 
is not accounted for in the A^(Hi) and Wco maps. 

In section 4.2.1 we have used the 7-ray emis- 
sivities per Hi atom and Wco unit to calibrate 
the CO-to-H2 conversion factor, following a well- 
established method. We can use a similar pro- 
cedure to correlate the 7-ray emissivities per Hi 
atom and per E(B — V)^.^,^ unit in the well resolved 
Gould-Belt clouds (see Fig. 4) where the spatial 
association between the Hi, CO and E(B — V)^.^^ 
maps allows locating the dark gas in the absence 
of kinematical information. 

The gHi. 1 and ^ebv emissivities found in the 
five energy bands exhibit a tight correlation 
(Fig. 13). As we did for CO, we fitted a linear 
relation, qebv = Xebv • 'Zhi, 1 + 9' using a max- 
imum likelihood method taking into account the 
errors and covariances of the emissivities. The re- 
sults are Xebv = (28 ±2) x 10^" cm^^ mag-^ and 
q = (-0.6 ± 0.2) X 10-^ cm-2 sr^^ mag-^ 
The good linear correlation implies similar spec- 
tra for the 7-ray emission from gas seen in the H i 
emission line and that associated with the excess 
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Table 3 

Masses for specific clouds, complexes or regions obtained from CO intensities and the 
Xco VALUES IN Table 2. For selected clouds we also report virial masses and, in the 
Gould Belt, the dark-gas mass obtained from the Xebv conversion factor determined 
IN 4.2.2. All masses are in units of 10^ Mq and the errors include only the statistical 

uncertainties on Xco or Xebv- 





/ 


/; 


(1 (kpc) 




il/vn('M) 


-^/v,r((/')) 




Cepheus 


[100, 117] 


[6,22] 


0.3"* 


0.37±0.02 


0.687 


0.903 


0.160±0.011 


Polaris 


[117,129] 


[18, 30] 


0.25'^ 


0.052±0.003 


0.208 


0.159 


0.031±0.002 


Cassiopeia 


[117,145] 


[2, 18] 


0.3^^ 


0.61±0.03 


0.893 


1.062 


0.34±0.02 


Gould Belt 


[100, 145] 


[-15,30] 


0.3 


1.47±0.08 








NGC 7538 


[107,115] 


[-5,5] 


2.65^= 


20±2 








NGC 281 


[120, 125] 


[-9,-5] 


3.0^ 


0.79±0.08 


1.205 


1.047 




Perseus arm 


[100, 145] 


[-10,10] 


3.0 


57±6 









^Grenier et al. (1989) 
''Heithausen & Thaddeus (1990) 
'^Moscadelli et al. (2009) 
<^Sato et al. (2007) 



Table 4 

TS = 2A(ln^) FOR the inclusion of the E(B - V)j.gg map in the fit in the different energy 

bands. 



energy range (GeV) 


TS 


0.2-0.4 


53.8 


0.4-0.6 


124 


0.6-1 


74.6 


1-2 


91.8 


2-10 


38.2 
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Fig. 13. — Correlation in the Gould Belt between 
the H I emissivities and the emissivities per unit of 
E(B — V)j.^g. Error bars show the statistical errors 
obtained on both emissivities in the five energy 
bands. The (red) line shows the best linear fit. 

reddening, thus confirming the need for normal 
additional gas to explain the correlated excess of 
both 7 rays and dust at the interface between the 
Hi and CO emitting phases of the Gould Belt 
clouds. 

Using the Xebv factor in these clouds, we can 
calculate the additional gas mass and compare it 
to the molecular mass seen in CO. We restrict the 
comparison to the CO mass, because the more dif- 
fuse Hi clouds are difficult to separate from the 
background H i disc extending to intermediate lat- 
itudes. To estimate the dark mass, we use only 
the positive residuals in the E(B — V)^.^,^ map. As 
discussed in section 2.3, the small negative resid- 
uals associated with the CO cores are likely re- 
lated with local variations in the dust tempera- 
ture or dust-to-gas ratio. The results are given in 
Table 3. The errors include only the statistical un- 
certainties on Xebv- The additional mass in the 
Gould-Belt clouds appears to be 40% to 60% of 
the CO-bright mass. We note that the sum of the 
dark and CO mass is closer to the virial one. We 
also note that EZRAS and SIMBA dust spectra in 
the Cepheus flare led to an independent estimate 
of its total mass, M = (0.43 ±0.18) x IO^Mq (Bot 
et al. 2007), which relates well with the total (CO 
plus dark) mass M = (0.53 ± 0.02) x IO^Mq we 
have obtained in 7 rays. 



5. Summary 

We have analysed the interstellar 7-ray emis- 
sion observed by the Fermi LAT in the region of 
Cassiopeia and Cepheus, successfully modeling the 
7-ray data as a linear combination of contributions 
arising from different gas complexes towards the 
outer Galaxy. 

The separation has allowed us to verify that the 
7-ray emissivity of local atomic gas is consistent 
with production by interactions with CRs with the 
same spectra as those measured near the Earth, 
but confirms the higher pion-decay contribution 
relative to some of the estimates in the literature, 
as found in Abdo et al. (2009c). This can be 
plausibly attributed to uncertainties in the local 
CR spectra, either in the measurement or from 
differences between the direct measurements and 
local interstellar space. 

Thanks to the correlation between an excess of 
dust and of 7-ray emission, with a spectrum equiv- 
alent to that found for the atomic and molecular 
gas, we have verified the presence of an excess of 
gas not properly traced by the standard Af(Hi) 
and Wco maps. In the nearby Gould-Belt clouds, 
the dark gas forms a layer between the H i and CO 
phases and it represents about 50% of the mass 
traced in the CO-bright molecular cores. 

The CR-density gradient in the outer Galaxy 
appears to be fiatter than expectations based on 
the assumption that CRs are accelerated by SNRs 
as traced by pulsars. It is also possible that the 
CR spectrum in the Perseus arm is harder than 
in the local arm. This hardening, which needs 
confirmation at high resolution with more LAT 
data to limit the potential contamination by hard 
unresolved point sources, could be linked to CR 
diffusion not far from their sources. 

We have measured Xco in several regions from 
the Gould Belt to the Perseus arm. The 7-ray es- 
timates are independent of the chemical and ther- 
modynamical state of the gas and also from as- 
sumptions on the virial equilibrium of the clouds. 
They correspond to a significant but moderate in- 
crease of Xqo with Galactocentric radius outside 
the solar circle, from (0.87 ±0.05) x 10^" cm^^ (K 
km s-i)-i in the Gould Belt to (1.9 ± 0.2) x 10^° 
cm~^ (K km s~^)~^ in the Perseus arm. 
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